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We establish the influence of nonlinear electron dynamics in the magnet ost at ic field of a hybrid 
planar free-electron maser on its gain and interaction efficiency. Even for the 'ideal' undulator 
magnetic field the presence of uniform longitudinal (guide) magnetic field potentially leads to the 
existence of chaotic zone around certain (magnetoresonant) value of the guide magnetic field. The 
width of the chaotic zone is given by the Chirikov resonance-overlap criterion applied to the normal 
undulator and cyclotron frequencies with respect to the coupling induced by the undulator magnetic 
field. Using analytical asymptotically exact solutions for trajectories of individual test electrons, we 
show that the magnetoresonant multiplier in electron trajectories is also present in the expression 
for the gain. The same Chirikov resonance-overlap criterion allows us to estimate analytically the 
maximal magnetoresonant gain of a hybrid planar free-electron maser showing that, in spite of 
the well-known drop in the gain for the exact magnetoresonance, the operation regime in the zone 
of regular dynamics slightly above the magnetoresonant value of the guide magnetic field is the 
preferable one. 

PACS numbers: 41.60.Cr, 05.45.-a 

Keywords: FEM/ubitron amplifier, chaotic dynamics and intcgrability, magnetoresonant growth rate 



I. INTRODUCTION 

Starting with the successful experiment [l], Q , previ- 
ous three decades have witnessed a spectacular devel- 
opment of theory and experiment of free-electron maser 
with guide (uniform longitudinal) magnetic field (hy- 
brid FEM). Utilizing the Dopplcr frequency upshift, free- 
electron masers and lasers posses a unique property to 
amplify and generate coherent electromagnetic radiation 
across nearly complete electromagnetic spectrum: from 
radio waves to vacuum ultraviolet [H . An interest to the 
hybrid FEM is reinforced by a promise to attain a sub- 
stantial microwave power level (up to a few gigawatts) 
in their planar configuration because of the use of su- 
perwide sheet electron beams (e.g., in Q emission and 
transport of a 140 cm-wide, 20 kA and 2 McV sheet 
electron beam for purported FEM applications was re- 
ported). Sheet electron beams allow one to weaken re- 
strictions on the maximal aggregate beam current posed 
by space charge effects while reaching the current values 
of about 30 -7- 50 kA. This makes the hybrid planar FEM 
one of the most attractive sources of powerful electro- 
magnetic radiation in the Terahertz Gap: the frequency 
range from 0.3 to 3 THz. 

In pioneering works 

@, B 0, H, HI analytical investi- 
gations of stationary regimes of microwave amplification 
and generation in a hybrid FEM were carried out; later 
these results were refined mainly through numerical simu- 
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lations. Experimentalists (lfj, [ll|, [l2j, [l3| reported a con- 
siderable loss of electron beam current and microwave 
power for a hybrid FEM for a certain range of values 
of the guide magnetic field. A study of chaotic parti- 
cle dynamics in free-electron lasers was undertaken in 
papers 0, [l5[, where the authors investigated effects 
of high-current (high-density) regime and the transverse 
spatial gradients in the applied wiggler magnetostatic 
field. Nevertheless, to the present day there still exist 
some fundamental questions, which either have no an- 
swers or only partial ones (cf. [H, pp. 430-439]): un- 
der what conditions motion of an individual test elec- 
tron becomes chaotic in the presence of guide magnetic 
field; what does influence the width of dynamical chaos 
zone around the magnetoresonant value (the undulator 
frequency is approximately equal to the cyclotron fre- 
quency) of the guide magnetic field; which maximal por- 
tion of initial kinetic energy from the longitudinal motion 
can be transferred to the transversally- vibrational motion 
(in an FEM the transfer of beam energy from longitudi- 
nal constant motion to the microwave field takes place 
indirectly through coupling of transversally-vibrational 
degrees of freedom to the microwave); does such an in- 
teraction with the microwave cause any widening of the 
chaotic dynamics zone; and, finally, what docs define the 
maximal value of the gain under the conditions of the 
magnetoresonance? 

We focus our attention on the operation of FEM 
around the magnetoresonant regime caused by the guide 
magnetic field usually used in FEM setups to enhance 
the efficiency of beam-microwave interaction and provide 
transverse confinement of electron beam. The study is 
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based upon analytical asymptotically exact expressions 
for electron trajectories developed recently [H[ as 
well as on the use of direct numerical simulations. More 
specifically, in the adopted approach we highlight the fun- 
damental role of nonlinear dynamical system describing 
motion of electrons in the combined magnetostatic spa- 
tially periodic (undulator) and uniform guide magnetic 
field. 

The next Section contains nonlinear and linearized in 
the microwave field self-consistent systems of equations 
governing a hybrid planar FEM amplifier. In Section HTT1 
we find asymptotically exact solutions of equations of mo- 
tion of an individual test electron in the magnetostatic 
field of a hybrid planar FEM and study their proper- 
ties. Section IIVI contains derivation of dispersion equa- 
tions and calculation of magnetoresonant growth rates. 
A comparison with direct numerical simulations of self- 
consistent nonlinear system of governing equations is pre- 
sented in Section [Vj The article ends with the Summary 
and Discussions Section. 



II. GOVERNING EQUATIONS 

We regard that at the entrance of the interaction re- 
gion (cross-section z = 0) of a hybrid planar FEM the 
electron beam is continuous and unmodulated. In the 
amplifier regime there also is a seed microwave signal. It 
can be shown that for a regular electrodynamics struc- 
ture the microwave field and electron beam current den- 
sity in the interaction region are periodic functions of 
time. Neglecting temporal harmonics generation, we will 
assume that there occurs resonant interaction only with 
one spatial harmonic of the microwave field. Under such 
assumptions the microwave field in an ideally conducting 
regular waveguide with the electron beam 



^(r,t)«Re{CS°(r,t)}=Re{C(z)e(f_ L )e- 4(wt - fc " z) }, 
B(r,t)«Re{CB°(r,i)}=Re{C(z)6(r L )e- 4(w '- fc " 2) } 



(1) 



differs from that of without the electron beam only by an 
amplitude multiplier C{z) (E°(f, t) and B°(f, t) obey ho- 
mogeneous Maxwell's equation); e(rj_) and b(fj_) are the 
membrane eigenfunctions of the coordinates f± = (x,y) 
in the waveguide cross-section S (microwaves and the 
electron beam propagate in the positive direction of the 
z-axis, see Fig. [l|; k z is the 'cold' propagation constant 
of the eigenwave. 

In the case of low-density electron beam, one can ne- 
glect its space-charge and, by averaging on the electrons' 
entrance phase, reduce the problem of interaction of elec- 
tron beam with a seed synchronous microwave of fre- 
quency to to the task of solving a single-particle equations 
of motion for arbitrary entrance time t e and Maxwell's 
equations in the form of Kisunko-Vainshtein's equations 
of excitation for regular waveguides H(| • It is exactly 
in this manner a highly successful th eory of the gyrotron 
was initially developed by Gaponov [2l| . 



Following [20, p. 200] (see also [22|), let us write the 
self-consistent nonlinear system of equations in the form: 



1 dt 

dC 
dz 



eE+-[vx (B p + B)}, 



(2a) 



7T fu 



(j ■ (E°f)e luJt dtdS, (2b) 



S—tt/u 



where 



j(r, t) = eV' 5(f— r k )i 



P o = g^Re J s (Wx) x b*(r ± )} ■ e z )dS. 

Here m e and e are the electron rest mass and charge, 
respectively; c is the speed of light; v = df/dt, 7 = 
(1 — v 2 / c 2 )^ 1 / 2 is the relativistic factor; j(r, t) is the elec- 
tron current density; r k (t) = r(t 7 t k ) and v k {t) = v(t,t k ) 
are the coordinates and velocity of a fc's electron, which 
entered the interaction region at the moment of time 
t*; the star (*) denotes complex conjugation; e z is the 



unit vector along the z-axis; the overbar in equation (|2b[) 
stands for the time average, which is effectively reduced 
to the averaging over the electrons' entrance phase. Ini- 
tial conditions for system @ consist in the specification 
of coordinates, velocities and the amplitude of microwave 
field at the entrance to the interaction region: r (t e , t e ) = 
(x,y,0), v(t e ,t e ) = (V x ,V y ,V l{ ) and C(0) = C . 

Magnetostatic field B p (r) consists of a guide (uni- 
form longitudinal) magnetic field B\\ and planar spa- 
tially periodic undulator magnetic field B w (r) = 
-B_L(0,cosh(27r[y - b' /2}/ \ w ) sm{2nz/ \ w ),smh(2n[y - 
b'/2]/X w ) cos(2ir z / \ w )) . Having the primary goal in 
demonstration of underlying physics, we intentionally 
consider a simplest possible model of a hybrid planar 
FEM (which accounts for the undulator through only one 
component of its spatially periodic magnetic field taken 
in the limit 2ir[y - b'/2]/\ w < I): 



(0,-B ± sin^z/AJ,-^), 



(3) 




FIG. 1: Schematic drawing of rectangular waveguide with 
mounted permanent magnet planar undulator immersed in 
the uniform longitudinal (guide) magnetic field. 
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where B^ is a constant amplitude [12|, [23[ and b' is the 
undulator gap width. However, the developed approach 
can be extended straightforwardly on 'realisable' undula- 
tor/wiggler magnetic fields and space charge dominated 
beams. 

In order to explore analytically the linear stationary 
regime of amplification, we linearize equations in the 
microwave fields E and B (B^ , B ± » c\E\/\v\, \B\). Tra- 
jectory of an electron is primarily given by its motion in 
the magnetostatic field 

r(t, t e ) ee f(t, t e ) = f Q (t) + Re{ri (*>**"<} (|f | » 

t = t — t e stands for the transit time of an electron. We 
also hold that C{z) is a slow function of z. In the zeroth- 
order in the microwave field approximation system ([2|) 
reduces to equations of motion in the magnetostatic field 



dt 



- [v x Bp(f )], 



(4) 



where v = dr Q /dt and 70 = (1 — Vq /c 2 )^ 1 / 2 . It then 
immediately follows that in this approximation the to- 
tal kinetic energy is conserved (i.e. 70 = const; Vq = 
v 2 (t e ,t e ) = V 2 + V 2 + V 2 ), which greatly simplifies the 
analysis of system (j4|). It should be emphasized that 
properties of solutions to equations ^ are fundamental 
to the consideration of microwave amplification and gen- 
eration in a FEM @, EH • In the first order (linear) 
approximation in the microwave field, we obtain a cou- 
pled system of equations for r\ and linearized equations 
of excitation 



m e 7o-~ 
dt 



Vl 



2^ (vo- Vl) 
To«o — -3 — 



zCF 



(5a) 



'-{[vi x B p (fo)] + Wo x (fi • V)B p (r)| r - =r - ]}, 



dC 
dt 



4P n 



(5b) 



where 



f(ro±,v Q ) = e(r Q j_) + c~ 1 [v xb(r a± )}. 



Here vi = dfi/dt, Io is the electron beam current at 
the entrance to the interaction region and F(t) is the 
microwave force estimated at the unperturbed electron 
trajectory. 



III. MAGNETOSTATIC RESONANCE: ENERGY 
TRANSFER AND DYNAMICAL CHAOS 



In this Section we examine dynamics of an individual 
test electron in the magnetostatic field © given by sys- 
tem ((4]). In components the equations have the form: 



i - u\\yo= - oj±z Q am(uj z /V\i), 
jjo + LJ\\X = 0, 
z = uj±xq sm(w zo/V\\), 



(6) 



where ujq = 27rV\\/\ w is the undulator frequency and 
wjii = \e\B± i|/(m e C7o); the overdots denote differen- 
tiation with respect to the transit time t. One can easily 
find two first integrals of this nonlinear dynamical system 
just by integrating its first two equations and take the en- 
ergy conservation law as the third linearly independent 



xq — u>\\yo — (2it) 1 uj±X w cos(ujoZo/V\i) = Fi = const, 
i/o + ojiiXq ee i<2 = const, Zq + Xq + j/q ee F a = const. 

Although, system ([6]) cannot be integrated in quadra- 
tures everywhere because it possesses non-zero Lyapunov 
exponents (see below Fig. [7] and [H, p. 208]). It should 
be noted that out of six quantities X w , B±, Bn, Vji, 
V x and V y only four their dimensionless combinations 
e = lo±/lo ee c/C/(7oV||) (K, is the conventional undu- 
lator parameter), <7 = u^/too, V x /V\\ and V y /V\\ define 
the behavior of nonlinear dynamic system Solution 
to system ( Qp c an be obtained by the use of a method of 
Lindshtedt |27l.l28j of asymptotic expansion of trajectory 
and frequencies in the small parameter e (see [l7| for a de- 



tailed exposition of the method in the case V x = V y = 0) . 
To the order o(e 2 ) the velocity components read: 



x (t) = Vx cos(Qy£- i>) + V\\ 



Vo(i) = -Vj_sin(ft||t-^) +«|| 



'0 

71q — fit 



cos(^ot) — COs(0||f) 



n 



sin(riiit) — — - sin(rioi) 
f2 



(7) 



zo(t) 



Vu — V1LO1 



-cos([f2o — tt\\]t + ip) cos([^o + Q\\]t — ip) 



2(fi -O||) 



2(fio + ^|| 



where V± = {V 2 + V 2 ) 1 ' 2 , sin^ = V y /V±, v\\ = kV\\, o(e 3 ) are given by the following self-consistent system of 
ilo = kujq and f2|| = trwy/co- Here k and a to the order 
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algebraic equations: 



Zk 2 + u 2 V; 2 (8 K 4 -13 K 2 a 2 +a 4 ) 



^ 2 (2K 4 + ll K 2 ( 7 2 -a 4 



V^n 2 



A(k 2 -<t 2 ¥ 



8V, 2 K 2 («; 2 - <7 2 ) 3 



8U 2 k 2 (k 2 



o- 2 ) 3 



a 



£ 2 (^ 2 + a 2 
4(k 2 - a 2 )'- 



-; (8) 



k and cr are found as functions of £ and <j$ with V x /V\\ and 
V y /V\\ as parameters. An analytical solution to Eqs. ([5]) 
can be found by successive iterations starting with 'non- 
renormalized' values 1 and <jq for k and cr, respectively. It 
is worth noting that Eqs. ([8]) order by order in e provide 
cancelation of secular terms in the solution procedure for 
system© (cf. [Tt}, which contains the limit of ©, Q 
and (|A.1[) for the case V x = V y = 0). For complete- 
ness, expressions for contributions to the velocities ([7]) 
of the order 0(e 2 ) are entered in the Appendix, Using 
Eqs. |(5J), and (|A.1[) one can check that Vq is con- 
served to the order o(e 3 ) and is equal to V 2 + V 2 + 
as required by the conservative nature of system ([6]). It 
should be also noted that this aggregate solution to the 
order o(e 3 ) is valid not only for the 'ideal' hybrid planar 
magnetostatic field ([3]) but also for its 'realizable' coun- 
terpart involving B w (r). The trajectories of electrons are 
then given by a straightforward integration. 

From ([7]) and (|A.1[) it follows that the motion of elec- 
tron is a superposition of constant motion with velocity 
uii and three-dimensional oscillations with normal undu- 
lator fio and cyclotron fi|| frequencies, which differ from 
their partial analogues luq and wu by 'rcnormalization' 
multipliers k and ct/cto- Recall that each subsystem of 
a nonlinear dynamical system possesses its own (partial) 
oscillation frequency if interaction between subsystems 
vanishes. Non-zero interaction modifies those oscillation 
frequencies by some amount [29l . |30| (represented here 
by multipliers k and a/ao); k also shows what portion of 
electron initial kinetic energy is transferred to the con- 
stant motion. 

As tends to fio the amplitudes of transversal and 
longitudinal oscillations increase. This signifies the ex- 
istence of an internal resonance (3l| . i.e. a resonant 
transfer of kinetic energy from longitudinal constant mo- 
tion to transversal and longitudinal oscillatory degrees 
of freedom. We call this situation magnetostatic reso- 
nance or magnetoresonance as is customary in the FEM 
theory [Hj]. 

A. Zero initial transversal velocity 

In the case when initial transversal velocity of electrons 
vanishes (V x = 1^ = 0), the expressions (?]), JSJ) and 
(|A.1[) are greatly simplified. Then, as mentioned above, 
the single-particle electron dynamics is defined only by 
two dimensionlcss quantities e and ctq. To calculate the 



normal undulator and cyclotron frequencies fio & n d fill, 
we need to solve for k and a the system of two algebraic 
equations ([5]) . To achieve maximal accuracy it turns out 
advantageous to solve this system numerically and use 
thus obtained values of k and a in the analytical calcu- 
lations. In Fig.[^]the results of such a solution are shown 
in firm lines. For comparison, the result of calculation 
of k and a via Fourier analysis of direct numerical solu- 
tion of Eqs. ([6]) for the velocity components is also given 
(shown by dots in Fig. In Fig. [3] the graphs of k 
and a vs. Bu are shown in more details for the guide 
magnetic field values around the magnetoresonance. At 
the magnetoresonance -B|| 6S « 6.67 kG the normal fre- 
quencies fio and fin experience jumps in order to comply 
with the conservation of kinetic energy in the dynam- 
ical system ([6]). As £ grows above a certain threshold 
value two zones of regular dynamics become separated 
by a chaotic zone centered about the magnetoresonant 
value of the guide magnetic field (see Fig. [5]). Typical 
Fourier spectra of velocity are shown in Fig. [4] Away 
from the magnetoresonant values of the guide magnetic 
field one can easily identify two distinct frequencies cor- 
responding to fio and fiy, their beating frequencies and 
harmonics. Around the magnetoresonance the spectrum 
is nearly continuous, which is one of the criteria of dy- 
namical chaos. 

Mutual orientation of the undulator and guide mag- 
netic fields causes the kinetic energy transfer from the 
longitudinal constant motion to undulator vibrational 
(frequencies luq and fio) and then to cyclotron vibra- 
tional (frequencies wy and fin) degrees of freedom of 
the nonlinear dynamical system ©; the coupling is pro- 
vided by the amplitude of undulator magnetic field B± 
(e = uj±/ljq). Such a picture allows one two possible ways 
of energy transfer: at magnetoresonant (fio ~ fiy) and 
non-magnetoresonant (fio 76 fi||) coupling regimes. Obvi- 
ously, the magnetoresonant regime requires smaller cou- 
pling and may become more advantageous for microwave 
FEM amplifiers and oscillators. It is then necessary to 
clarify how accessible are these regimes of FEM opera- 
tion. 

Viewing dynamical system Q as a nonlinear pendu- 
lum with external modulation of its vibration amplitude 
(cf. the last equation in ([6])), one can estimate the condi- 
tion of appearing of chaotic layer near the separatrix. An 
analytical estimate is consistent with numerical calcula- 
tions and coincides with the Chirikov resonance-overlap 
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FIG. 2: k and cr as functions of the guide magnetic field Bu . 
K and and a are depicted in black and red colors, respectively; 
for reference (To(-Bii) is given by the chain line. The darkened 
area corresponds to the zone of chaotic dynamics as given by 
the Chirikov resonance-overlap criterium (|9}. 



FIG. 3: K and a as functions of the guide magnetic field for 
-By values around the magnetoresonance. k and and a are 
depicted in black and red colors, respectively; for reference 
cro(Bii) is given by the chain line. 



criterium [33| : 

|fio-fi||| <u± (|«-tr| <e), (9) 

i.e. appearance of the chaotic state takes place whenever 
the difference between normal frequencies of the system 
becomes less than the coupling. As will be shown in Sec- 
tion [IV] the absolute value of the ratio ui±/(£lo — fin) 
(called below the magnetoresonant multiplier) defines 
the microwave gain. Results of numerical calculations 
of its maximal value for zones of regular dynamics and 
practically accessible e(B±,V\\) are given in Fig. [5] (for 
each value one needs to optimize B\\ ; solid and dotted 
lines correspond to approach of the chaotic zone from 
low and high values of the guide magnetic field, respec- 
tively). This shows that in the magnetoresonant case 
there exists limitation on the fraction of energy trans- 
ferred from constant motion to oscillatory (transversal 
and longitudinal) degrees of freedom posing a fundamen- 
tal limit on an FEM efficiency. A calculation shows that 
Xo(£) + Vo(t )/zq(1) w 0.48 (bars denote time averaging). 
In Fig. [S] we present a contour plot of k = v\\/V\\ as func- 
tion of dimcnsionless coupling e and guide magnetic field 
do- The region of chaotic dynamics is given in black. In 
Fig. [7] a map of the major Lyapunov exponent as func- 
tion of s and tro is shown. The calculation included a 
check that the sum of all Lyapunov exponents is equal to 
zero as mandatory for a conservative dynamical system 
because of the Liouville theorem. 



B. Non-zero initial transversal velocity; 
chaos control 

Let us start with demonstration of some features of so- 
lutions ([7]) for a non-zero initial transversal velocity. As 
follows from (J7J) and ([5]) , one can have either suppression 



(V x > 0) or enhancement (V x < 0) of the velocity ampli- 
tude of cyclotron vibrations, respectively. In particular, 
for V y = and V x = v\\uj±Qo/(QI — fijj) there occurs a 
complete, i.e. to the first order in the small parameter s, 
suppression of cyclotron vibrations in the transverse to 
the z-axis plane. In Fig. [8] this situation is demonstrated 
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0.5 1 1.5 2 2.5 
Normalized Frequency 0)1 0)^ 

FIG. 4: Logarithmically scaled Fourier spectra of xo as func- 
tion of normalized frequency uj /coo ■ Top figure shows regular 
dynamics while bottom figure corresponds to chaotic one. 
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FIG. 5: Maximal value of magnetoresonant multiplier as func- 
tion of dimensionless coupling e. 



for V y = and both signs of V x of the same absolute 
value \V X \ = 0.1c, V\\ = 0.66c, B x = 1.11 kG, the value 
of Bu is chosen to be 6.67 kG and corresponds to the 
magnetoresonance for V x < (only projections of elec- 
tron trajectories in the xy-plane are shown). 

It should be noted that in general the presence of non- 
zero initial transversal velocity influences the shape of 
chaotic dynamics zone. The influence of V y is symmetric 
regarding its sign and, hence, less relevant while that 
of V x is much more important (see ([8])), In particular, 
the x-componcnt of initial velocity of a certain sign can 
either suppress (V x > 0) or enhance (V x < 0) the zone of 
chaotic dynamics as shown in Figs. [9] and 1101 Non-zero 
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FIG. 6: Contour plot of k = «||/^| as function of dimension- 
less coupling e and guide magnetic field <ro- The region of 
chaotic dynamics is given in black. 



V x also changes position of the edge of regular dynamics 
zone at zero guide magnetic field. This is due to the 
fact that under such conditions position of separatrix on 
the phase plane zq-zq depends on the value and sign of 
V x and is independent of V y . Position of the edge at 
zero guide magnetic field (oo = 0) is calculated to be 
e=([l + (V x /V ll T} 1 /* + V x /V ll )/2. 

Having considered separately the particular cases of 
zero and non-zero initial transversal velocity, below we 
will study some generic properties of nonlinear dynamical 
system ([5]). First, it should be emphasized that nonlinear 
dynamical system ([5]) is degenerate in the sense that at 
zero amplitude of undulator magnetic field or zero guide 
magnetic field the dynamics is characterized by only one 
normal frequency instead of two. Under such conditions 
(after subtraction of constant motion along the z-axis), 
electron trajectories constitute invariant curves, while in 
the generic case they wind up (surfaces diffeomorphic 
to) two-dimensional invariant tori (see below). Second, 
an important feature of solutions in the cases of zero and 
non-zero initial transversal velocity consists in cancela- 
tion of non-physical secular terms in the series expres- 
sions. This is due to the use of a method of Lindshtcdt 
for finding electron trajectories. Third, as mentioned at 
the beginning of this Section, nonlinear dynamical system 
© possesses three linearly independent first integrals Fj, 
(i = 1, 2, 3) that are not all in involution 

{F 1 ,F 2 } = -u h {F 1 ,F 3 } = {F 2 ,F 3 } = 0. (10) 

({Fi,Fj} denote the Poisson brackets). Therefore, it 
does not fall under conditions of the Liouville theorem. 
However, observing that the identity, /, is also a neces- 
sary member of the Lie algebra of first integrals of system 
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FIG. 7: Map of the major Lyapunov exponent. 
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FIG. 8: Projections of electron trajectories in zy-plane. (a) 
V x — —0.1c; (b) V x = 0.1c. Bounding box is of a standard 
waveguide size a x b = 8.6 mm x 4.3 mm. 



© by virtue of the first relation in (JTDJ) , one can 
infer that the condition of non-commutative integrabil- 
ity holds for the system under investigation (for details 
see [HI, P- 190]). It then follows that we can apply the 
Nekhoroshev theorem to find the dimensions of invariant 
tori (the loci of the electrons' trajectories), which foliate 
the system's phase space, [Hj|. There exist four first 
integrals (I, Fi,i — 1,2,3) for a three dimensional 
nonlinear dynamical system ([H]), two of which (/ and F%) 
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FIG. 9: Contour plot of k e Vy/Vy as function of e and ao'- 
suppression of zone of chaotic dynamics (V x — 0.3 Vj| ) . The 
region of chaotic dynamics is given in black. 



constitute a complete set of independent first integrals 
in involution. Then, the dimensions of invariant tori 
are equal to two. Under natural projection into the 
system's configuration space one can visualize such 
two-dimensional tori. Namely, in the reference frame, 
which moves with the mean velocity of electron (see 
Eqs. O and (|A.1|> ). 7y = (1 - u 2 /c 2 ) -1 / 2 , its trajectory 
winds up a (surface diffcomorphic to) two-dimensional 
torus (Fig. fTT|) . Such a regular behavior is exhibited 
by the dynamical system under investigation only for 
particular range of dimensionless parameters e, uq, 
V x /Vii and Vy/V|| (cf. [36|) as can be understood from 
Figs. O ED and [T0J When parameters change in such 
a way that the motion of electrons becomes chaotic, 
we can observe the disintegration of the corresponding 
invariant torus shown in Figs. [T!?l and [l~3l In Fig. [T3il the 
trajectory of electron is drawn for a small period of time 
and looks like winding up a two-dimensional torus, while 
for a longer period of time it does not (Fig. [iU]). 

Summarizing this Section, we observe that there are 
two regimes for pumping kinetic energy into transversal 
degrees of freedom of a hybrid planar FEM: the first one 
(advocated here as more, if not the only, efficient one 
for the terahertz waveband) consists in low amplitude of 
undulator magnetic field and the guide magnetic field as 
close to its magnetoresonant value as only consistent with 
the regular dynamics; the second regime (usually consid- 
ered in the literature) would also provide the same ratio 
LU±/(flo — through a high absolute value of |fio — 
given by an off-resonant guide magnetic field -By and a 
necessarily high value of amplitude of undulator magnetic 
field B±. It, however, seems that values of B± in excess of 
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FIG. 10: Contour plot of k = U||/V]| as function of e and ao: 
enhancement of zone of chaotic dynamics (V x = — 0.3Vj| ) . 
The region of chaotic dynamics is given in black. 
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FIG. 11: Electron trajectory (winding of two-dimensional 
torus) under conditions of regular dynamics. 




Y\\(za~V\\t) (cm) 0.5 <T = 0.77 

FIG. 12: Disintegration of invariant torus for a short period 
of time. 



10 kG arc currently practically unattainable while guide 
magnetic fields of up to 50 kG can be presently achieved. 



IV. LINEAR AMPLIFICATION: 
MAGNETOSTATIC RESONANCE AND 
MAXIMAL GAIN 

As mentioned in Sec. [D] analytical knowledge of elec- 
trons' trajectories ([7]) in magnctostatic field allows 
one to develop a gyrotron-like analytical linear theory 
of microwave amplification. In particular, in Ref. [37[ 
such an approach was constructed for interaction of an 
electron beam with a mode of an arbitrary regular waveg- 
uide in a hybrid FEM. It should be noted that analyt- 
ical theory of this type is valid not only far away from 
the magnetoresonance (^o ~ fill, cf. [38}) but also in its 
close vicinity. Here we will apply this theory to study 
amplification of microwaves in a hybrid planar FEM un- 
der conditions of magnetoresonance using as an exam- 
ple amplification of TE 10 and TE 01 modes of rectangular 
waveguide with wide a and narrow b sides. Exactly these 
modes allow one to utilize the advantage of enhanced 
current of sheet electron beams (cf., for example (fT4|) ). 
Specifically, in each case we need to solve system of 
ordinary differential equation with quasi-periodic coeffi- 
cients. Its exact solution is achieved by application of a 
quasi-periodic analogue of Floquet's theorem and, under 
some additional conditions, comes in the form of infinite 




FIG. 13: Disintegration of invariant torus for a long period of 
time. 



series in undulator and cyclotron harmonics, [39l . pp. 90- 
96]. Under particular synchronism conditions the series 
may be truncated to retain solution pertinent to this 
very synchronism neglecting other combination ones (see 
[20I l2lj] ) . In this case, we can write C(z) = Co exp (iSk z z) 
for the microwave amplitude (\Sk z \ <C \k z \), where 5k z is 
the spatial growth rate. 

A. Undulator synchronism 

Let the undulator synchronism condition, [l8l ]. is ful- 
filled 

u - k° z v\\ w O . (11) 
I 



ri 



^001 



Zl 



and Su> = uj — k z vu — Qq is a small mismatch from the 
ideal synchronism (Slo = 0). It then follows that a ma- 
jor contribution to the spatial growth rate of the waveg- 
uide mode exists even for an electron beam with zero ini- 
tial transverse velocity. For conciseness we will assume 
V±_ = down to the end of this subsection. On the one 
hand for V± = expressions (fl2|) describe a well-known 
result (see [4(|) that interaction between the transverse 
component of oscillatory motion and microwave results 
primarily in axial bunching, which becomes the main 
source of instability (if V± equals to zero then radial 



Using the fundamental TE10 mode (k® = [(ui/c) 2 — 
(tt/o) 2 ] 1 ^ 2 ) of a rectangular waveguide as an explicit 
demonstration example, we solve (|5a| to the order o(e 2 ) 
(retaining only summands up to the second order in 
V_l/v«) and obtain 



n(t) = rie r + r o 0ie0 + z x e z , (12) 



where 



and azimuthal bunching is weaker because in this case 
ri and ro4>i are inversely proportional to the small quan- 
tity (5k z v\\ — Stu) while z x is inversely proportional to 
its square). On the other hand (unlike in [4(|) we in- 
fer this from analytical expressions (|12p . which are valid 
for all regular dynamics parameters as found in Sec. IIIII 
and, specifically, near the magnetoresonance. Substitut- 
ing (TT2"]) in (|5bp . one then gets the dispersion equation 
for the TE10 mode. 

Results of the outlined calculation in the case V x = 
V v = are as follows: 



rA^lelCo lu^ eicp(i[6k z V\\ — Su>]t) j . Ini Z \ 11 \ - , V v 



< sin(fi||i) cos 



■r 



4m e c7o fig — Q 2 Sk z v^ — 8ui 1 Vet 

27r 2 c 2 (n - \a [ + n l{ J J 16a 2 n 2 v 2 ° [ 11 ^ C ° S \a [ OyVJ' 

iir\ w \e\C u±Sl exp(i[5k z Vu - Su]i) r / w r _ , V y 

— < cos(S2iiij cos I — \x 1 



cos(f2||i) cos^— [a 



4m e c7o fig — Q 2 5k z v\\ — 5uo I \a Q 

a\ w Sl {l\\(u> + tt\\)V± . /7r r _ U yl \ A^Q2y2 _ ^ 



iTra\e\ujCo Wj_f2|| cxp(i[<5fcjV|| — <$u;]t) j A^f^u; + 0|j)Vj \ . f n r- Vy 
2m e c^ n 2 -Q 2 (Sk z v l{ - Suj) 2 I 1 16a 2 r2jj7j[^2 J sm Va^ + 



fill 



- ^) 2 

V|| 



w 1 



2cC7 7o 3 (^-fi 2 ) 2 



fi 2 w 



■ sin 



! (^), TE 10 mode = [(c/c) 2 - (n/a) 2 } 1 / 2 ); 



abc- Ka_J' (13) 
^sin 2 (^), TE 01 mode = [(c/c) 2 - (tt/6) 2 ] 1 / 2 ); 



where U = m e v 2 /(2\e\) is the non-relativistic beam volt- two complex conjugated roots of each of Eqs. (flU)) con- 
age of constant motion. The conditions for existence of 
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sist in positivity of their respective cubic discriminants. 
These yield in each case the thresholds for such instabil- 
ities to occur, pill ]. 

For the TEio mode such a threshold condition reads 



5u> 



> 



3 

27o 



1/3 ( ^11 A 2/3 -2/3 (™\ 

\abc 2 UJ \aJ- 



For the ideal synchronism (Suj = 0) the spatial growth 
rate of the TEio mode is given by 



Im<5fc, 



_ V3 / tuj|/o| y/ 3 / ^±»ii \ 2 / 3 

~2~7o" \2abc 2 UJ 



O 2 



fl 2 



sin 



2/3 



(?)■ (») 



Re <5fc 2 provides a correction to the cold propagation con- 
stant, k®, caused by the presence of electron beam. Re- 
sult (TTj)) is similar to those found previously (cf., e.g. (I^ l 
and also implies that there could exist a substantial en- 
hancement in the (microwave) gain, G = 8.63Im<5fc 2 , 
because of the presence of guide magnetic field. It is 
also significant that analytical expression (TTJJ provides 
a close approximation of the growth rate, which is in a 
good agreement with the direct numerical simulations of 
the nonlinear self-consistent system of relativistic equa- 
tions of motion and equations of excitation ©. This is 
mainly due to the fact that in J7o and Q\\ we account for 
the initial electron velocity and magnitudes of magneto- 
static fields not only through the definitions of luq and u>\\ 
but also via 'rcnormalization' multipliers k and a (cf. © 
and Fig. [2]). More importantly, utilizing results obtained 
in Sec. Mil we can provide an analytical estimate for the 
maximal attainable gain. To accomplish this task let us 
substitute from © to find 



O 2 



O 2 

HI 



< 



f2n + Oil 



(15) 



Recall that there exists two ways of having the value of 
magnetoresonant multiplier, luj_/(Qo — fin), close to 1: 
one should either work far from the magnetoresonancc 
(and provide for a strong undulator amplitude B±) or, 
as advocated at the end of Sec. IIII1 operate near the 
magnetoresonancc f2o ~ ^|| with a moderate undula- 
tor amplitude (see also Fig. [S]) . From the inspection of 
the right hand side of (JTSJ) , we observe that its value is 
about 1/2 around the magnetoresonance. This provides 
another reason for preferring a hybrid planar FEM opera- 
tion regime slightly above the magnetoresonance in order 
to ensure the regular dynamics of electron motion as de- 
tailed in Sec. IIIII (see also Figs. [2J [3] and [5]). We, therefore, 
find for the upper limit on the maximal (near) magne- 
toresonant gain of a hybrid planar FEM under conditions 
of ideal undulator synchronism (Scu = 0) and Vj_ = 



C res ^3.7/HW/ 8 c 

max 7o ^ abc 2 U J 



sin 



2/3 



(?)• <»> 



The same conditions of regular dynamics lead to an im- 
portant result that around the magnetoresonance the 



maximal gain is independent of B± (although entering 
Vu in U through n, it cancels out completely by virtue of 
synchronism condition = lo / (k° z {uj) + 2n/X w )); as a 
function of frequency, G™| x - (k° + 2tt/ Xyj) 2 / 3 ^ 1 ^- 1 / 3 . 

Under the assumption of undulator synchronism in the 
limit -B|| = (a = and k 1 — 3e 2 /4) an electron beam 
without initial transverse velocity does not interact with 
the TE m0 modes of a rectangular waveguide. Then the 
TEoi mode turns out to be the lowest one amplified by 
such an electron beam, and the bottom line of Eqs. (fT3"|) 
can be used to estimate analytically the gain of a planar 
FEM amplifier without the guide magnetic field (cf. Q) 



G\ 



B||=0 



Although, it should be noted that magnetoresonant gain 
(£lo « and ui±/(flo — fin) 1) calculated using the 
bottom line of Eqs. (fl"3"|) is always greater than that one 
given by Eq. (fT?)) . 

Another important characteristic of an amplifier is its 
tunability. In the case of undulator synchronism the fre- 
quency tuning of an FEM amplifier is achieved by chang- 
ing initial axial velocity, Vji , of the electron beam and 
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FIG. 14: In the top figure results of analytical calculations of 
the maximal resonant gain (solid line) and initial axial veloc- 
ity of the beam (dashed line) for the fundamental TEio mode 
are presented. In the bottom figure we show the gain (solid 
line) and initial axial beam velocity (dashed line) at Bu = 
for the TEoi mode. All quantities are shown as functions of 
amplified frequency. 
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turns out to be limited only by the requirement of single- 
mode operation regime. In Fig. [T?] for the fundamental 
TEio mode one can see that the maximal resonant gain 
depends weakly on the frequency in the operating range. 
To provide the single-mode operation regime and interac- 
tion only with the forward wave, the frequencies are cho- 
sen to range from 18 GHz to 32 GHz in this case. Results 
of gain calculations for interaction with the TEoi mode 
are also shown in Fig. [Mj The initial axial velocity must 
also change with the frequency (v\\ = u/ (k z (cu) + 2n j 'X w )) 
as shown in the both parts of Fig.[l4]to maintain the ideal 
undulator synchronism (in the top figure for calculation 
of Vj| from a given vu we choose K = 0.8, see Fig. [6] in 
the bottom figure to obtain V\\ from a given v\\ one needs 
to solve the equation v\\ = k(V\\ )Vji , see (jHJ and note that 
in this case a = 0, V x = V y = 0). 



B. Hybrid synchronism 

Under the condition of hybrid (undulator-cyclotron) 
synchronism, fl8j . 

w-k°v\\ wn +fi||, (18) 

we again can solve equation (|5aj) to the order o(e 2 ) re- 
taining only summands up to the second order in Vl/um 
and those which contain the square of the small quan- 
tity [8k z v\\ — 5uj) in the denominators. The result for the 
fundamental TEio mode has again form (|12|) but now 
with 



0, ro0i = - 



\l\e\n*c 0l *vl 



cj ± exp(i{[5k z v\\ — Sw]i — ip}) 



32 



7rm e c 3 7o(l + v \l\ /c 2 ) v f\ — fl\\ {$kzV\\ — Suj) 2 



3a 



k° z X wK vpfin 



n\ w \e\QoCoV± 



K + a 27t(k + ct) 4c 2 ( 
ex.p(i{[Sk z v\\ — 5u>]t — ip}) 



r \cos[ - a 

k — a) j \a 



8m e C7o(l + vp 2 /c 2 )v\\ fl - 



{- 



-9 1 



2k a — 4k<t 



(5k z v» — Stu) 2 

2 '2a 3 2 7 2 (^ 3 



6n 2 a + 8na 2 - a 3 ) 



27r(K + cr) 4c 2 L k(k - <t){k + 2a) 



7 2 k(k + a)(n + 2a) 



} cos G 



(19) 



YlL 

fin 



and7|| = (1 — V 2 /c 2 )~ 1 / 2 . It then follows that the bunch- 
ing mechanism is substantially axially-azimuthal, which 
is in contrast with the undulator synchronism where it 
is predominantly axial. In the non-rclativistic limit the 
azimuthal mechanism does not contribute to the bunch- 
ing being a manifestation of dependence of cyclotron fre- 



quency on the Lorentz factor (cf. [2l|). Substituting the 
obtained result in (fT2")) and then in (|5b[) . one again gets 
the dispersion equation for TE 10 mode under the hybrid 
synchronism condition. 

In this manner, we obtain 



Sk z (Sk z - 6 -^) 2 
v \\ 



8cu l0 (i + vfif/c 2 )(n 2 -n 2 ) 2 



irVfuj 
a 3 bc 

ab^il 2 



Qi 



'(" 

v a 



Q 



TE . 



[y 



fin 



, TEio mode; 
, TEqi mode; 



(20) 



which are expansions to the second order in V_\_/v\\ of the Here the analytically exact expressions for the relativistic 
bulky analytical totally relativistic dispersion equations. factors read 
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-)TE 

no 



= 1 



V\l\ C $7T(T 



7 2 (k + ct)(k 3 



_ » (k - g) 
2lok 

- 2k 2 o - 4kct 2 



2cj 
2d 3 ) 



272(K-Cr)K 2 (K + 2(T) 



-jTE 
^01 



= 1 - 



2irv 2 j 2 a 27rwj|7Q6cr 



(ft + Cr)(K 2 + K<7 — (T 2 ) 

k 2 (k + 2a) 
-I 7Tu|7^ (k + cr)(J 

27r 2 cf|746cr 3 



Note that fc° = [{uj/c) 2 - (tt/o) 2 ] 1 / 2 for the TEi mode 
and fc° = [(w/c) 2 — (tt/^) 2 ] 1 ^ 2 for the TEoi mode, respec- 
tively. 

In the non-relativistic limit (substituting = 1 

and keeping only the top line of expression (|20[) with 
70 = 7|i = 1 and v\\/c <C 1) the obtained dispersion equa- 
tion for the fundamental TEio mode is similar to that one 
known in the literature (321 . Eq. (18)]. However, even in 
this limit, not only the use of Eq. ([20)) is justified by 
our procedures for all parameter values compatible with 
the regular dynamics (i.e. near the magnetoresonance) 
but it also provides a good quantitative agreement with 
the direct numerical simulations of self-consistent non- 
linear system ^ because of the 'renormalization' mul- 
tipliers k and a (e.g. if one takes k to be equal to 1 
then 1 + ?; 2 7 2 /c 2 = r ) 2 ). It also features analytically cal- 
culated dependence of the right hand side of dispersion 
equation on the entrance position x of electrons (the fac- 
tor cos 2 (7r[x + V y /£l\\]/a)). Thus, one can observe that 
there exist two equivalent optimal entrance position of 
electrons, e.g. for iV^I/f^i < a the planes x = — V y /Qu 
(for V y < 0) and x = a — V y /Q\\ (for V y > 0), placed 
near the maxima of gradient of the microwave electric 
field (since the electrons should not touch the waveguide 
walls the optimal entrance positions are not always pos- 
sible). This, in principle, allows one to double the elec- 
tron current by using two electron beams (with the same 
x-componcnts and equal in value but opposite in direc- 
tion y-components of initial velocity) and enhance the 
total spatial growth rate by a quarter (while simultane- 
ously providing for reduction of aggregate space-charge 
effects). 

Near the magnetoresonance, analogously to the case of 
undulator synchronism, there occurs a substantial reduc- 
tion in the mean velocity of constant motion, vu , because 
of kinetic energy transfer to mainly transversal oscillatory 
degrees of freedom. This leads to the gain growth, which 
is only bounded by transition to the chaotic state (see 
Fig.rrO BT, es = 3.65 kG). As seen in Fig. [17] the magne- 
toresonance also influences essentially the dependence of 
ideal synchronism (Stu) frequency on the guide magnetic 
field. It should be also emphasized that, as we found in 
Sec. IIII B[ the sign of V x has a substantial impact on the 
location of zones of regular and chaotic dynamics, there- 
fore, this sign can also influence the microwave amplifi- 
cation. Unfortunately (because the problem involves at 



least six parameters) , to find a set of paprameters which 
maximizes the right hand side of (|20p . we are bound to 
confine themselves to detailed numerical calculations. It 
turns out that the maximal spatial growth rate in the 
regular dynamics zone is attained near the magnetores- 
onance for V x > and V y > 0, i.e. suppression of dy- 
namical chaos near the magnetoresonance (achieved by 
a positive V x , cf. Fig. [5]) may allow one to enhance the 
spatial growth rate. 

For completeness, we also provide the result of calcu- 
lations for the TEoi mode given by the bottom line of 
Eq. (|20p . although to implement such an interaction in 
a single-mode regime one needs to take additional efforts 
for an efficient mode selection (see, e.g., [![)• 



V. NONLINEAR SIMULATIONS OF 
MICROWAVE AMPLIFICATION 

For numerical simulations of nonlinear regime of am- 
plification it is advantageous to rewrite equations of mo- 
tion (|2a[) taking as independent variables the coordinate z 
and time of entrance of electrons to the interaction re- 
gion t e (the time of arrival t = t(z, t e ) of electrons, which 
entered the- interaction region at the time t e , to the cross- 
section z becomes a dependent variable) in the form 

dz V m e c7 / dz 

dt m e r y dr± p± 
dz p z dz p z ' 

where E(r±,z,t), B(r^,z,t) and B p (r±,z) are given 
by (TTJ) and ([3]), respectively; f± = (x,y); it is also more 
convenient for numerical calculations to introduce the 
momentum of an electron p = m e v'y instead of its ve- 
locity (pj_ = (p x ,p y ), 7 = [1 + p 2 /(m 2 c 2 )} 1/2 ). The ini- 
tial conditions then are x(z=0,t e ) = x, y(z=0,t e ) = y, 
t(z=Q,t e ) = t e and p(z=0,i e ) = m e {V x , V y , Vj|) 7 o. Simi- 
larly, using the charge conservation law and the fact that 
in the stationary regime electrons, which enter the in- 
teraction region at the time t e separated by an integral 
multiple of the period of the amplified microwave, go 
along identical trajectories, we can write the equations 
of excitation (|2b[) for a thin electron beam as follows [43|, 
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p. 31]: 



dC 

dz 



2tt 



4irP 



(22) 



Here the initial condition for the amplitude C(z) reads 
C(z=0) = Co. In the discrete model of electron beam we 
assume that over the period of the amplified microwave 
particles enter to the interaction region in regular time 
intervals 2w/(u}N), hence, the numerical finding of so- 
lution to nonlinear self-consistent system (|21[) and (j2"2"]) 
consists in the simultaneous solving of 6iV + 1 first-order 
nonlinear ordinary differential equations. 

Numerical simulations of FEM characteristics are usu- 
ally concerned with those, as a rule integral, quantities, 
which can be directly measured experimentally (e.g. sat- 
uration power, efficiency, start current, etc.). However, 
a substantial advantage of numerical calculations also 
lies in the opportunity of detailed reconstruction of elec- 
trons dynamics and their interaction with the microwave 
field in an FEM. Thus, according to analytical results 
of Sec. IIV Al the predominant phasing mechanism on the 
undulator mechanism is the axial one. A numerical verifi- 
cation to this fact is provided in the top of Fig. [15] In the 
region 1, out of the uniform at entrance {z = 0) electron 
beam, there occurs bunch creation under the influence 
of the seed microwave in its decelerating phase. The in- 
teraction of electron beam with the microwave is almost 
linear and the microwave power P(z) = \C(z)\ 2 Pq grows 
exponentially with 2lm5k z = 0.024 cm -1 . In the region 2 
the electron beam is maximally bunched and efficiently 
amplifies the microwave (simultaneously about 0.35N of 
the beam electrons, which entered the interaction region 
at z = in the accelerating phase of microwave, draw the 
power from the microwave); the total microwave power in 
this region increases almost linearly. In the region 3 the 
electrons, which entered the interaction region at z = in 
the decelerating phase of microwave, do not interact with 
it, but the microwave power growth takes place because 
of interaction with approximately 0.5./V of beam electrons 
that entered the interaction region at z = in the acceler- 
ating phase. In the region 4 the energy flow from the elec- 
trons, which amplify the microwave, and those that draw 
microwave power reaches the balance and the microwave 
power attains saturation. In the region 5 there occurs 
the second (like the region 2) consolidation of bunches 
at z = 230 cm, which subsequently leads to the second 
maximum of microwave power at z = 280 cm and so on. 
It is also worth noting that electrons, which entered the 
interaction region at z = in the accelerating phase of 
the microwave, are weaker trapped than those entering 
the interaction region at z = in the decelerating phase. 
The small-amplitude fast vibrations shown in Fig. 1 151 are 
caused by the oscillations of longitudinal velocity of elec- 
trons around its mean value v\\ (see j7j and (|A.1|) ). Note 
that parameter values used to produce Fig. [TS] are typi- 
cal for experimental designs under development, but it is 
not optimal in terms of efficiency of microwave amplifica- 
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FIG. 15: Results of numerical calculations of electron beam 
bunching and microwave power amplification. 



tion because of weak effective pumping of electron oscil- 
lations. As found in Sec. IIV1 the optimal amplification is 
achieved near the magnetoresonance fio ~ ^11 • This op- 
erational regime of a hybrid planar FEM is studied in the 
literature analytically and numerically in comparatively 
less detail, mainly, inasmuch as the major attention of 
researches has been devoted to hybrid helix and coaxial 
FEM schemes. For example, in the hybrid helix scheme, 
one uses an annular electron beam, which adiabatic en- 
trance is incompatible with the magnetoresonant condi- 
tion [i^l since the electrons entering the interaction re- 
gion with different radial separations from the symmetry 
axis of the undulator magnetic field reach orbits charac- 
terized by a substantial mean velocity spread (this orbits 
fail to be the desired stationary hclicoidal ones) and do 
not provide an adequate amplification of microwaves. In 
the planar FEM on the undulator synchronism this draw- 
back is offset substantially because a sheet electron beam 
enters in the symmetry plane of the undulator magnetic 
field and, therefore, can efficiently amplify microwaves 
under the magnetoresonant condition. 

Results of numerical calculations of various character- 
istics of hybrid planar FEM as functions of the guide 
magnetic field Bu are plotted in Figs. [16] and [171 One 
can see that the optimal amplification frequency (the fre- 
quency providing for the maximal growth rate) both for 
undulator and hybrid synchronisms (see (jlip and (|18p ) 
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FIG. 16: In the top figure solid black (gain G) and red (fre- 
quency /) lines are obtained from analytical expressions (|11[) 
and (|13[) : the dots of respective colors are the results of direct 
numerical simulation of Eqs. (|21[) and (122[) . In the bottom 
figure black (efficiency r/) and red (saturation length L aat ) 
dots are obtained through numerical simulation of Eqs. (|21[) 
and (|22l) . All quantities are calculated as functions of Bit 
for the TEio mode amplification under the undulator syn- 
chronism. The darkened area corresponds to the analytically 
calculated zones of chaotic dynamics around the magnetores- 
onance [B\ ea « 6.67 kG). 



FIG. 17: In the top figure solid black (gain G) and red (fre- 
quency /) lines are obtained from analytical expressions (|18p 
and (|20|) : the dots of respective colors are the results of direct 
numerical simulation of Eqs. (|21[) and (|22[) . In the bottom 
figure black (efficiency rf) and red (saturation length L sa t) 
dots are obtained through numerical simulation of Eqs. (|21[) 
and (|22|) . All quantities are calculated as functions of B« 
for the TEio mode amplification under the hybrid synchro- 
nism. The chain lines show the position of magnetoresonance 
B^ es m 3.65 kG (analytical calculations indicate an exponen- 
tially thin zone of chaotic dynamics for the chosen parame- 
ters). 



depends on the guide magnetic field since the mean ve- 
locity of constant motion. vu , is also a function of Bu at 
least through the 'renormalization' multiplier n. Under 
the undulator synchronism the gain, G, and efficiency, 
77, attain their maximal values (G ~ 0.47 dB/cm and 
?7 w 5.8%) for the guide magnetic field (-By ~ 7.5 kG) 
greater than its magnetoresonant value. A theoretical 
estimate of the width of chaotic dynamics zone for the un- 
dulator synchronism gives B\\ £ (6.6, 7.3) kG around the 
magnetoresonance B^ es ps 6.67 kG. The situation is simi- 
lar under the hybrid synchronism with a slight difference. 
The gain has a pronounced maximum G ~ 0.37 dB/cm 
for the guide magnetic fields greater than but close to 
the magnetoresonant value B^ es ss 3.65 kG, however, 
the maximal efficiency 77 us 1.75% is attained at a small 
finite separation (Bu w 4.9 kG) from the magnetoreso- 
nance. Both cases of undulator and hybrid synchronisms 
show a larger efficiency of microwave amplification for 
guide magnetic fields greater than their respective mag- 



netoresonant values since under such conditions a larger 
portion of the energy of constant longitudinal motion is 
transferred to the transversal oscillations of electrons. 

It should be also noted that whereas the growth rate 
is a slow function of the amplitude of undulator mag- 
netic field B± (or, equivalently, of e), the efficiency 77 
grows almost linearly with it up to a certain critical 
value r\ cr ~ 20% (e.g., e cr ~ 0.2 under the undula- 
tor synchronism). If e exceeds this critical value then 
the coherent amplification of microwave signal in the 
(near) magnetoresonant regime becomes impossible. In 
Sec. IIIII we established that because of the Nckhoroshev 
theorem the motion of electrons in the magnetostatic 
field of hybrid planar FEM is characterized by invariant 
tori, overwhelming majority of which (according to the 
Kolmogorov-Arnold-Moser theory, see |30() under small 
perturbations of motion by a microwave signal will not 
be destroyed but will only be slightly deformed. Thus, 
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interaction of electron beam with the microwave signal 
does not lead to a notable change of zones of chaotic and 
regular dynamics if this interaction is sufficiently small 
(cf., e.g., Figs. 151 and 1161 for the undulator synchronism). 
However, if the interaction becomes larger of a certain 
critical value in the efficiency (a measure of beam-wave 
interaction) and, therefore, apparently, larger of a certain 
critical value of e then such a notable change of zones of 
chaotic and regular dynamics does take place, i.e. there 
occurs a strong widening of the domain of disintegrated 
invariant tori. 



VI. SUMMARY AND DISCUSSIONS 

We applied Kisunko-Vainshtein's equations of excita- 
tion for regular waveguides to the development of self- 
consistent analytical linear theory of a hybrid planar 
FEM amplifier, which is valid not only far away but 
also around of the magnetoresonant value of the guide 
magnetic field (f2|| f^o). Nonlinear numerical simula- 
tions were undertaken to clarify the validity of the linear 
approximation. Although, as mentioned at the begin- 
ning of this article, a number of approaches to the linear 
theory of a hybrid FEM has been already formulated, 
a new one presented here allowed us both to provide a 
consistent and unified analytical description of electron 
dynamics for all values of the guide magnetic field and, 
more importantly, to find an efficient characterization of 
the dynamical chaos present in the system. In particu- 
lar, conditions of suppression of the dynamical chaos are 
formulated and the influence of hybrid planar FEM op- 
erational parameters on these conditions is clarified. It is 
also worth mentioning that the zones of suppressed beam 
transport found in [12j | are in good correspondence with 
our (semi-)analytical calculations for the zones of dynam- 
ical chaos (cf. Fig. 1 in [l^] and Figs. 121 and 1101 of this 
article). This is achieved because the method of Lind- 
shtedt employed by us to obtain test individual electron 
trajectories in the magnetostatic field of a hybrid planar 
FEM is capable of high precision in analytical calcula- 
tions. It should be also mentioned that this method is 
sufficiently general to be applicable to describe motion of 
charged particles in spatially inhomogeneous static mag- 
netic fields like, for example, that in jl2l. [23j . 

The origin of dynamical chaos is connected to the 
possibility of onset of stochastic layer around the sep- 
aratrix in the system phase space (under the condition 
uj ± /lu = [(1 + V^/V 2 ) 1 ' 2 + Vx/V^/2 motion of an elec- 
tron in the undulator magnetic field takes place along 
the separatrix) caused by the presence of guide magnetic 
field, which plays the role of non-trivial perturbation ca- 
pable of destruction of this separatrix (see Figs. 
and UOp . The magnitude and sign of the x-component 
of initial velocity of electrons influence the position of 
separatrix in the phase space and allows one to adjust 
location of zones of regular and chaotic dynamics when 
the guide magnetic field is present. An approximate con- 



dition that determines the onset of chaos is given ana- 
lytically and found to be in a good quantitative agree- 
ment with numerical simulations. We interpret this con- 
dition as the Chirikov resonance-overlap criterion, i.e. 
the chaotic behavior occurs whenever the absolute value 
of the difference between the normal undulator, and 
normal cyclotron, f2||, frequencies becomes less than the 
coupling, luj_, induced by the undulator magnetic field. It 
seems also that the analytical approach developed in this 
paper has a strong potential for a development of analyt- 
ical nonlinear in the microwave signal three-dimensional 
theory of a hybrid FEM. 

The transfer of kinetic energy of an electron motion 
between longitudinal and transversal degrees of free- 
dom is studied and it is found that the maximal frac- 
tion, (fy± — l)/(7 — 1), of the mean transversal kinetic 
energy is less than 30% for the case of regular dynamics. 
Taking into account that transversal degrees of freedom 
of electrons are the source of energy for the microwave 
field one can claim that the maximal efficiency of a hy- 
brid planar FEM may not exceed these 30%. Using the 
Nehoroshev theorem we showed that motion of an in- 
dividual test electron in the magnetostatic field is in- 
tegrable and characterized by two-dimensional invariant 
tori for some range of parameters of hybrid planar FEM. 
As a result from the Kolmogorov-Arnold-Moser theory it 
then follows that for relatively moderate values of the mi- 
crowave and space-charge field the electron trajectories 
will stay regular and amplification can still be accom- 
plished. Through numerical simulations, we determined 
that the broadening of chaotic region for electron motion 
does not occur if the interaction between the microwave 
and electron beam is not too large (e.g., if such a measure 
of interaction intensity as the efficiency is less than 20% 
on the undulator synchronism). 

From linearized in the microwave field equations of mo- 
tion and excitation dispersion equations are derived for 
the undulator and hybrid synchronisms. We showed ana- 
lytically that around the magnetoresonance on the undu- 
lator synchronism the gain is nearly completely indepen- 
dent of the amplitude of the undulator magnetic field. 
This circumstance is known in the literature but only 
as a result of numerical simulations (cf. Fig. 2 in [HI). 
Physical origin of this effect lies in the fact that the gain 
is a function of amplitude of transversal oscillations of 
electrons in the magnetostatic field of a hybrid planar 
FEM. Around the magnetoresonant value of guide mag- 
netic field this amplitude, according to the chaotization 
criterium ([5]), turns out to be bounded and is indepen- 
dent of the amplitude of undulator magnetic field. On 
the undulator synchronism the obtained analytical ex- 
pression for the gain provides values close to the results 
of nonlinear numerical simulations. Such an accuracy is 
achieved because of the high-precision analytical calcula- 
tion of the unperturbed by the microwave proper frequen- 
cies and electron trajectories. For the hybrid synchro- 
nism we also obtained the dispersion equation but in this 
case the analytical expression for the gain in the zone of 
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regular electron dynamics immediately below and above 
the magnetoresonant value of the guide magnetic field 
does not provide values close to the results of numerical 
calculations. We attribute this discrepancy to a greater 
amplitude of oscillations of the longitudinal electron ve- 
locity and, therefore, to a greater deviation (comparing 
to the case of the undulator synchronism) of the resultant 
system of ordinary differential equations from that with 
the strictly periodic coefficients. Solutions to such or- 
dinary differential equations with the quasi-periodic co- 
efficients does not follow patterns of their counterparts 
with the periodic coefficients especially in the vicinity of 
(combination) resonances between their proper frequen- 
cies. Nevertheless, we showed analytically and verified 
through numerical simulations that operation of a hybrid 
planar FEM in the (near) magnetoresonant regime is op- 
timal in order to achieve the maximal gain and efficiency. 
The interaction between transversal degrees of freedom of 
electrons and microwave field under such conditions turns 
out to be maximal out of all other possible operation pa- 
rameters of a hybrid planar FEM on each of the synchro- 
nisms (experimentally for a hybrid planar FEM oscillator 
such an operational regime was established, for example, 
in [Hj])- The major result here is that one needs not only 
to have a small value of e in order to maintain the rela- 
tion v\\ 3> v± (v± is the mean transversal velocity) but 
also to hold the ratio e/\k — <r\ = oj±/\flo — fiii| as close 
to the unity as possible thus providing for the maximal 
gain and efficiency. 

Here we have not considered any specific arrangements 
usually applied for a smooth (at best adiabatic) entrance 
of an electron beam to the interaction region, which lead 
us to neglect of the velocity and position spread of elec- 
trons at the entrance to the interaction region. The 



developed techniques are fully capable of treating these 
effects but taking them into account in the framework 
of the current consideration would make our treatment 
overcomplicated and hide behind the technicalities the 
physics underlying the nature of dynamical chaos in a hy- 
brid planar FEM. This paper also does not deal with the 
influence of the space-charge field on operational charac- 
teristics of a hybrid planar FEM. It is necessary to notice 
that space-charge field strongly decreases the efficiency of 
interaction between an electron beam and amplified mi- 
crowave if the regime of operation is positioned far away 
from the magnetoresonance. However, there are strong 
indications that not only the amplification of a microwave 
signal by an electron beam becomes maximal around the 
magnetoresonance but also the defocusing influence of 
the present space-charge field turns out to be minimal. 
A detailed investigation of the influence of potential (ir- 
rotational) and rotational parts of the space-charge field 
on the operation of a weakly-relativistic hybrid FEM will 
be published elsewhere. 
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APPENDIX: VELOCITY TO THE ORDER o(e 3 ) 

Here we present terms of the order 0(e 2 ), which should 
be added componentwise to expressions (JT)) to obtain so- 
lutions to the electron velocities valid to the order o(e 3 ). 
These additional contributions read 
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